A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

03(11.01.2; 11.17.3; 12.03.1; 13.18.1; 13.20.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 



An evolutionary model for GHz Peaked Spectrum Sources 
Predictions for high frequency surveys 

G. De Zotti 1 , G.L. Granato 1 , L. Silva 2 , D. Maino 2 , and L. Danese 2 

1 Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 1-35122 Padova, Italy 
, 2 SISSA, International School for Advanced Studies, Via Beirut 2-4, 1-34014 Trieste, Italy 
ON ' 

ON \ Received 19 July 1999 



o 



> 

(N 

oo 

(N 
(N 

ON 
ON 

Oh 
i 

o 



Abstract. We have explored, in the general framework of 
the "young source" scenario, evolutionary models for GHz 
Peaked Spectrum (GPS) galaxies and quasars which re- 
produce the observed counts, redshift and peak frequency 
distributions of currently available samples. Substantially 
different cosmological evolution properties are found for 
the two populations: the quasar luminosity function must 
evolve strongly up to z ~ 1, while the data on galaxies 
may be consistent with no evolution. The models show 
that GPS sources (mostly quasars) may comprise quite a 
significant fraction of bright (S > 1 Jy) radio sources at 
v > 30 GHz if the intrinsic distribution of peak frequen- 
cies extends up to ~ 1000 GHz. In any case, however, their 
fraction decreases rapidly with decreasing flux and their 
contribution to small scale fluctuations in the frequency 
range covered by the forthcoming space missions MAP 
and Planck Surveyor is expected to be minor. 

Key words: Galaxies: active - quasars: general - Ra- 
dio continuum: galaxies - cosmic microwave background 
- Submillimeter 



1. Introduction 

One of the major uncertainties affecting estimates of mi- 
crowave fluctuations due to discrete extragalactic sources 
(Toffolatti et al. 1998; Sokasian et al. 1999; Gawiser et al. 
1999) arises from the very poor knowledge of the abun- 
dance of radio sources with strongly rising spectra in the 
frequency range from a few tens to a few hundred GHz, 
where the most sensitive experiments to map primordial 
anisotropies of the Cosmic Microwave Background (CMB) 
are carried out. 

GHz Peaked Spectrum (GPS) radio sources [see O'Dea 
(1998) for a comprehensive review] selected at v < 5 GHz 
have a fairly flat distribution of peak frequencies, suggest- 
ing the existence of an hitherto unknown population of 
sources peaking at mm wavelengths (O'Dea & Baum 1997; 
Crawford et al. 1996; Lasenby 1996). Grainge & Edge 
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(1998) report the detection of 50 GPS sources brighter 
than 50 mJy at 5 GHz and spectra still rising above 
10 GHz. The observed peak frequency of several of these 
sources was found to be > 43 GHz; one has a redshift of 
3.398 so that the emission peak is above 190 GHz in the 
rest frame. 

It is thus possible that GPS sources significantly con- 
taminate experiments aimed at obtaining very high accu- 
racy (at the /iK level), high resolution maps of the CMB, 
such as those expected from NASA's MAP and ESA's 
Planck Surveyor missions. 

Unfortunately the very limited information currently 
available makes it very difficult to carry out reliable pre- 
dictions of the confusion noise due to these sources. On the 
other hand, physical models, although still schematic, may 
provide a useful guide for such predictions. An investiga- 
tion of this kind may also be useful in connection with the 
on-going searches of GPS sources peaking at high frequen- 
cies (Guerra et al. 1998; Grainge & Edge 1998; Cooray et 
al. 1998). 

Two main scenarios have been proposed to explain the 
properties of GPS sources. One hypothesis is that their 
ages are similar to those of classical double radio sources 
and are kept compact for a large fraction of their lifetimes 
by interactions with dense gas in their environment (van 
Breugel et al. 1984; O'Dea et al. 1991; Carvalho 1998). 
Alternatively, GPS sources may correspond to the early 
stages of evolution of powerful radio sources (Phillips & 
Mutel 1982; Carvalho 1985; Fanti et al. 1990, 1995; Read- 
head et al. 1996a,b; O'Dea & Baum 1997). The latter pos- 
sibility is currently favoured because the media surround- 
ing radio sources appear to be insufficiently dense to in- 
hibit the radio source growth for a large fraction of its 
lifetime (Fanti et al. 1995). Also, direct evidences of very 
short kinematic ages have been found by Owsianik & Con- 
way (1998), Owsianik et al. (1998, 1999) and Stanghellini 
(private communication) . 

In this paper we adopt the young, evolving sources 
scenario. We have elaborated (Sect. 2) on the analytic 
self-similar evolution model proposed by Begelman (1996, 
1999) to account for the main observed population prop- 
erties of GPS sources (Sect. 3). In Sect. 4 we present pre- 
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dictions for counts of GPS sources at MAP and Planck 
Surveyor frequencies. Our main conclusions are summa- 
rized and discussed in Sect. 5. 



2. Outline of the adopted evolution model 

We make the following assumptions: 

1. The initial radio luminosity function (in units of 
Mpc~ 3 dlogii -1 ) is described by a power law: 



-^z,min — -^i — ^i,maxi 



(i) 



where Li is the luminosity before absorption. 

2. In the GPS phase, the properties of the sources are de- 
termined by the interaction of a compact, jet-driven, 
overpressured, non thermal radio lobe with a dense 
interstellar medium (Bcgelman 1996, 1999; Bicknell 
et al. 1997); the timescale of the interaction is very 
short in comparison with the cosmological-expansion 
timescale, so that the luminosity evolution of indi- 
vidual sources occurs at constant z. As the radio 
lobe expands in the surrounding medium, the emit- 
ted radio power decreases with the source age, r, as 
Li oc r~ n , and its linear size I increases as / a r e . 
If the density of the surrounding medium scales with 
radius as p e oc r~ n , we have (Begelman 1996, 1999) 
r\ = (n+4)/[4(5— n)] and e = 3/(5— n). There is a clear 
anticorrelation between intrinsic turnover frequency, 
Up, and linear size (O'Dea & Baum 1997): v p oc l~ s , 
with 8 ~ 0.65. It follows that v v scales with time as 
v p oc t~ a , with A = Se. 

3. The spectra of GPS sources are described by: 

if v < v v 



L v 



f „ j (yl v v) ' " 

^ p X 1 {v/v p )- a ifv>v 1 



(2) 



v 



with a a = 0.8 and a = 0.75, the mean values found by 
Snellen et al. (1998b). 

As the radio lobe expands, the peak luminosity L p varies 
as the consequence of two factors: the decrease of the emit- 
ted radio power and the decrease of v p . Hence: 



Lpivp) — L p .iT p — Lp 



- P /\ 



= L r 



7)1 A — a 



(Vpjvp .i) _ ' 1+1 ^ A ' ) , the epoch dependent luminosity func- 
tion at a given frequency v (Mpc -3 d log L^ 1 GHz -1 ) 
writes: 



n{L u ,Vp,z) = n (l + z) 3 x 

Lp^i{L u , Vp) 
L*{z) 



with Lp.i(z) oc Li(z) and p = —r\ + a\. 

If the birth rate of GPS sources is constant on 
time scales much shorter than the cosmological-expansion 
timescale, the peak luminosity function per unit dlogi p 
is: 



n(Lp) oc lVp. 



(4) 
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(—) 




\Vp,i) 



-(1+1/A) 
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where L±(z) is the redshift-dependent normalization lu- 
minosity. We have assumed luminosity evolution and 
adopted a very simple parametrization for it: 



' (1 + z c ) k if z> z c 



(6) 



The redshift z c at which luminosity evolution levels off is 
a model parameter. We have normalized monochromatic 
luminosities to L = 10 32 ergs -1 Hz -1 . 

The luminosity function at a frequency v 2 is related to 
that at a frequency v\ by: 

!{v 2 /v l y aa if vi < v 2 < v p 

(vi/v p ) a "-(v 2 /vp) a if vi < v p < v 2 (7) 
(v 2 /vi) a if v p < vi < v 2 

The number counts per steradian of GPS sources brighter 
than S v at the frequency v, with an observed peak fre- 
quency max(i/, ^ Pimin ) < v Pi o < ^ p ,max are given by 

r min[z f ,z m (S u )] 

N (> S v ] Vp.a > v) = I dz — x 



o 



dz 



x 



min[i/ p>max (l+«),^ p> i] 



dv p x 



'max[i'p ;m i n (l+2),i/(l+2)j 
/■log L max (z) 

x / d log L v n(L v ,u p ,z), 

J log L mia {S„,z,v p ) 



(8) 



where Zf is the redshift of formation of the first GPS 
sources (we have set Zf — 3.5, the maximum observed 
redshift of a GPS source), z m is the maximum redshift at 
which sources can have a flux > S v , i m ; n is the minimum 
luminosity of a source of given z and v p yielding a flux 
> S v , dV/dz is the volume element within a solid angle w 
in a Friedman- Robertson- Walker universe (A = 0): 



(3) 



dV _ c 

~dz~ ~ lh u> JT+ zf{\ + n z y/ 2 

2-0 



d L = —z 1 + z 



S v 



H a ~ \~ ' "2 + Oz + 2(l + fiz) 1 /2 
L v K{z) 



Also, since, in this case, the number of sources of age r 

within dr is simply proportional to dr and dr/dv p oc LC(z) = (1 + z)- 



(9) 

(10) 

(11) 
(12) 
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Similarly, the number counts of GPS sources with an ob- 
served peak frequency v p , m in < v Pi q < min(V, u p>max ) are 
given by 

rmin[z f ,z m (S v )] ^y 

N (> S„; v pfi <v) = J dz — x 

rmin[v (l+z),^ p ,max(l+z)^p,i] 

X / dv p X 

"'^P,mi„(l + Z) 
flogLmax(z) 

x / d\ogL v n{L v ,v p ,z) , (13) 

./log L m i n (S l/ fZjVp) 

Throughout this paper, we adopt H = 50 kms -1 Mpc -1 
and fi = 1. 

The distribution of observed peak frequencies per unit 
dv Py o in a flux limited sample, A/"(f Pj o; > S v ), is given by: 

rmin[z f ,z m (S v ,v Pi0 )\ ^y 

JV (v p ,o; > S„) = J dz — x 

/log imax(2^p,o(l+z)) 
d log L v n[L v , z, v pfi (l + z)](14) 
-Jgimin(S 1/ ,Z,l/p, (l+z)) 

3. Comparison with observations 

3.1. Samples 

The available information mostly comes from the sam- 
ples by Snellen et al. (1998b) and by Stanghellini ct al. 
(1998). Marecki et al. (1999) have recently defined a rela- 
tively large sample, but still with a limited fraction of op- 
tical identifications and redshifts. A search of GPS sources 
peaking above 10 GHz has been carried out by A.C. Edge 
and co-workers; however only a very preliminary report 
has been published so far (Grainge & Edge 1998). 

The sample by Snellen et al. (1998b) comprises 47 
sources selected at 325 MHz with a limiting flux den- 
sity of approximately 18 mJy; it was also required that 
,S(5GHz) > 20mJy. Whenever possible, an inverted spec- 
trum between 325 and 609 MHz was used as a criterion 
to select candidate GPS sources; otherwise, the selection 
was based on the 325-5000 MHz spectral index. The 325- 
609 MHz selection yielded 14 sources over an area of 119 
square degrees; this is probably a lower limit to the true 
areal density of GPS sources since the 5 GHz flux limit 
may introduce a bias against GPS sources with lower peak 
frequencies and steep spectra beyond the peak. The 325- 
5000 MHz selection yielded 33 sources over 522 square de- 
grees; the lower areal density, compared with the previous 
case, although not very statistically significant, may reflect 
a stronger bias against lower peak frequencies. The condi- 
tion that spectra of candidate GPS sources can be fitted, 
within the observed frequency range, by a self-absorbed 
synchrotron spectrum, discriminates against sources with 
high values of u pfi (> 8.4 GHz). Snellen et al. (1998a) 
obtained optical identifications and R and/or I magni- 
tudes for 41 of them. Redshifts were measured by Snellen 



et al. (1999) for 19 of these. Since GPS galaxies show a 
tight apparent magnitude-redshift correlation (Snellen et 
al. 1996), rcdshift estimates can be derived from photo- 
metric data; best fit relationships are provided by Snellen 
et al. (1996). 

Stanghellini et al. (1998) have defined a complete sam- 
ple of 33 GPS sources brighter than 1 Jy at 5 GHz, over 
an area of approximately 24600 square degrees. All objects 
are optically identified and have magnitude estimates; red- 
shifts are available for all but 5 of them. This sample is 
complementary to that of Snellen et al. (1998b) not only 
because of the selection at a much higher frequency and 
at brighter radio fluxes, but also because, with only four 
exceptions, the peak frequencies of sources are below the 
frequency of selection, while the opposite is true for the 
other sample. 

The sample by Marecki ct al. (1999) comprises 76 
sources with 5 GHz fluxes > 200 mJy over an area of about 
7655 square degrees. The GPS classification, however, may 
be doubtful for 28 of these sources whose spectral shape is 
not well determined. As mentioned by the authors them- 
selves, a further uncertainty arises from the variability, 
typical of flat spectrum sources, which may also lead to 
a misclassification of some sources as GPS's. Only 21 ob- 
jects are identified (3 galaxies and 18 quasars) and have 
measured redshift. 

The sample selected by Edge and co-workers is partic- 
ularly well suited to investigate the abundance of sources 
peaking at very high frequencies which could be tricky 
contaminants of CMB anisotropy maps. Grainge & Edge 
(1998) report the detection of 50 GPS sources brighter 
than 50 mJy at 5 GHz and peak frequencies above 10 GHz 
over an area of about 2000 square degrees; they point out 
that some GPS sources may have been missed because of 
various selection biasses and, therefore, their result should 
be regarded as a lower limit. No details on the redshift and 
peak frequency distributions are given. 

3.2. Constraints on model parameters 

GPS galaxies and quasars turn out to have different red- 
shift, rest-frame peak frequency, linear size and radio mor- 
phology distributions (Stanghellini et al. 1996; Snellen et 
al. 1998a; Stanghellini et al. 1998). Thus, they must be 
treated as different populations. Hence, we have fitted sep- 
arately the redshift and v p $ distributions of GPS galaxies 
and quasars in the samples by Stanghellini ct al. (1998) 
and Snellen et al. (1998b). The counts reported by Grainge 
& Edge (1998) were used as an additional constraint. The 
best fit values of the parameters were derived using the 
routine "amoeba" (Press et al. 1992) exploiting the down- 
hill simplex method in multidimensions. 

It became clear quite soon that there is no need for 
cosmological evolution to account for the observed red- 
shift distributions of GPS galaxies. Therefore, to minimize 
the number of free parameters, a no-evolution model was 
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Table 1. Normalization constants and model parameters 





Galaxies 


Quasars 


L (ergs" 1 Hz" 1 ) 




10 32 


u p>i (GHz) 




1000 


z f 




3.5 


no [Mpc" 3 (dlogL,,)" 1 GHz" 1 ] 


5.0 x 10"* 


6.1 x 10~ 13 


L min (300 MHz)/L 


6.7 


0.36 


Lmax(300 MHz)/L 


2.9 x 10 5 


4.2 x 10 s 





0.75 


0.59 


n 


10.7 


9.6 


A 


8.7 


3.9 


Zc 




1 


k 




12.2 



adopted = Lq in Eq. da)]. Only very few "local" 

GPS sources are known; therefore the local luminosity 
function cannot be determined directly. The power law 
luminosity function adopted here is fully characterized by 
four quantities: the normalization uq, the slope (5, the min- 
imum and maximum luminosities. Additional model pa- 
rameters are rj and A. 

Table [l] gives the best fit values of the parameters for 
both galaxies and quasars. The normalization of the lu- 
minosity function and the minimum and maximum lumi- 
nosities refer to the frequency of 300 MHz. 

Figs. and | compare the observed peak frequency 
and redshift distributions of GPS galaxies and quasars in 
the samples by Stanghellini et al. (1998) and Snellen et al. 
(1998b), respectively, with those yielded by the model. 

Redshifts of the 5 sources (four galaxies and 1 of uncer- 
tain identification, assumed to be a galaxy) in the sample 
by Stanghellini et al. (1998) lacking spectroscopic mea- 
surements were estimated using the magnitude-redshift 
relationships derived by Snellen et al. (1996). 

We have assumed that the 19 objects in the sample 
by Snellen et al. (1998b), identified as star-like by Snellen 
et al. (1998a), are quasars. As suggested by Snellen et al. 
(1998a), we have adopted the faint source, rather than the 
quasar, as the identification of B1647+6225. Of those, 12 
have measured redshifts (Snellen et al. 1999); their red- 
shift distribution was taken as representative of the full 
sample (in practice, the observed number of quasars in 
each redshift bin was multiplied by 19/12). Only 6 of the 
28 "galaxies" have a spectroscopic redshift; for the other 
alleged galaxies redshift estimates were derived from the 
magnitude-redshift relationships. 

The agreement between model and observed distribu- 
tions is very satisfactory: the formal value of \ 2 per degree 
of freedom is ~ 1 both for galaxies and for quasars; a word 
of caution is in order, however, in view of the many uncer- 
tainties following from the various selection effects men- 
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Fig. 1. Peak frequency and redshift distributions of 
quasars and galaxies in the sample by Stanghellini et al. 
(1998). The solid lines show the observed distributions. 
The dotted lines those implied by the model. 



tioned in the previous subsection (very difficult to quan- 
tify accurately) from photometric redshift estimates (es- 
pecially for the sample by Snellen et al. 1998b) and from 
some doubtful identifications. 

The model also yields 56 GPS sources (47 quasars 
and 9 galaxies) brighter than 50 mJy at 5 GHz and with 
Vpft > 10 GHz, over an area of 2000 square degrees, to 
be compared with a lower limit of 50 sources (no further 
specifications given) reported by Grainge & Edge (1998). 

Assuming that the range of values of isq (0.4 GHz < 
v < 10 GHz) listed in Table 4 of Marecki et al. (1999) cor- 
responds to the range of peak frequencies of GPS sources 
selected by these authors, the model yields 44 GPS sources 
(26 galaxies and 18 quasars) satisfying their selection cri- 
teria. As mentioned above, although the Marecki et al. 
(1999) sample comprises 76 sources, only 48 of them have 
an unambiguous GPS classification. 

We stated earlier that no evolution of GPS galaxies 
is required by the data. On the other hand, the derived 
evolutionary properties of GPS quasars are quite extreme, 
although the increase with redshift of the space density at 
any given luminosity is moderated by the flatness of the 
local luminosity function and by the relatively low value 
of the redshift at which the increase stops. 

The local luminosity functions of galaxies and quasars 
have almost identical slopes, significantly flatter than 
those of conventional radio sources (Dunlop & Peacock 
1990). The coincidence of the slopes of local luminosity 
functions of the two GPS populations may well be for- 
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Fig. 2. Peak frequency and redshift distributions of 
quasars and galaxies in the sample by Snellen et al. 
(1998b). 

tuitous. In the case of galaxies, the flat slope is required 
to reproduce redshift distributions extending to z > 3 
without resorting to cosmological evolution. In the case of 
quasars, a flat slope may be reminiscent of the effect of rel- 
ativistic beaming on luminosity functions (Urry & Shafer 
1984; Urry & Padovani 1991). Indications that quasar ra- 
dio flux densities are moderately increased by Doppler 
boosting (while GPS galaxies are not) are discussed by 
O'Dea (1998). 

An independent estimate of the local luminosity func- 
tion of GPS sources was recently reported by Snellen & 
Schilizzi (1999) who also found it to be flatter than that 
of steep spectrum radio sources. They argue that this is 
consistent with a scenario whereby GPS sources evolve to 
Compact Steep Spectrum (CSS) sources and to large-scale 
double radio sources, provided that the radio luminosity 
increases with time in the GPS phase and decreases in the 
subsequent phases. 

Since only qualitative arguments are given by Snellen 
& Schilizzi (1999) it is difficult to analyze quantitatively 
their model. We may mention, however, that we were un- 
able to obtain a consistent fit to the data sets considered 
here with parameter values implying an increase of the 
radio power with increasing age of GPS sources. On the 
other hand, according to Eqs. (2) and (3), the monochro- 
matic radio luminosity, L v , of GPS sources at frequencies 
v < v v varies as L v cx r v v~ aa cx t-p~ Aq ». Inserting the 
values of parameters given in Table 1, it turns out that, 
in the case of galaxies (but not in the case of quasars), L v 
does increase with time (cx t 2 8 ). 



Fig. 3. Predicted counts of GPS galaxies (left-hand 
column) and quasars (right-hand column) at three 
Planck Surveyor frequencies for different ranges of ob- 
served v p fi (GHz) and a maximum intrinsic value u Pl i = 
1000 GHz. 

The parameters r/ and A, characterizing the evolu- 
tion with the source age of the emitted power and of the 
turnover frequency, are much larger than expected in the 
framework of the self-similar evolution model by Begelman 
(1996, 1999). The physical meaning of this result needs to 
be further investigated. 

4. High frequency counts of GPS sources 

The high frequency counts of GPS sources further depend 
on the maximum value, v Pt i, of the rest-frame peak fre- 
quency. We have considered two extreme cases. The results 
shown in Fig. || for three Planck Surveyor frequencies cor- 
respond to a very large value, v v/l = 1000 GHz. On the 
other hand, since there is direct evidence that values of 
Up > 190 GHz do exist (Grainge & Edge 1998), we may 
take 200 GHz as a lower limit to v Pi i . The corresponding 
counts are shown in Fig. ^. 

The results confirm the rough preliminary estimate by 
De Zotti et al. (1998): from several tens to a few hundreds 
of such sources could be detectable by Planck Surveyor 
instruments. 

Not surprisingly, given the high luminosities of GPS 
sources at low frequencies and their inverted spectra, 
sources with the highest peak frequencies (mostly quasars, 
whose peak frequency distribution is weighted towars 
higher values, as shown by Figs. |l] and ||), show up at 
the brightest flux levels. If the rest-frame peak frequen- 
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Fig. 4. Predicted counts of GPS galaxies (left-hand 
column) and quasars (right-hand column) at three 
Planck Surveyorfrequencies for different ranges of ob- 
served Vpfi (GHz) and a maximum intrinsic value v V i = 
200 GHz. ' 



cies can reach very large values (up to ~ 1000 GHz), high- 
Vp t o sources are expected to dominate the bright end of 
GPS source counts up to ~ 100 GHz. A comparison of 
the results shown in Fig. || with the estimates by Toffo- 
latti et al. (1998) and Sokasian et al. (1999) suggests that 
GPS sources may in this case, correspond to a substan- 
tial fraction (30 to 70%) of bright radio sources (> 1 Jy) 
at v > 30 GHz. On the other hand, since their model 
counts are rather flat, the fraction of GPS sources rapidly 
decreases with decreasing fluxes; at 0.1 Jy their fraction 
should drop to ~ 5%, implying that they are a minor con- 
tributor to small scale fluctuations due to extragalactic 
sources. 

The areal density of sources peaking at very high fre- 
quencies obviously decreases rapidly with decreasing v v .i. 
For example, the 30 GHz counts of GPS quasars brighter 
than 1 Jy with observed v p .q > 30 GHz decreases by about 
a factor of 4 as u Pt i decreases from 1000 to 200 GHz 
(see Fig. [!]), while the corresponding total counts of GPS 
sources decreases by about a factor of 2. 

Given their rather low areal density, only very large 
area surveys would allow the selection of significant sam- 
ples of these sources. Planck Surveyor's all sky surveys are 
unique in this respect. 



5. Conclusions 

In the general framework of the "young source" scenario 
we have worked out exploratory models reproducing the 
observed counts, redshift and peak frequency distribu- 
tions of currently available samples of GPS galaxies and 
quasars. The derived values for parameters characterizing 
the evolution with the source age of the emitted power 
and of the turnover frequency, however, are substantially 
larger than those usually quoted. Therefore, the physical 
interpretation of these results needs to be further investi- 
gated. 

With the due allowance for the many selection effects 
and uncertainties related to the identification of bona fide 
GPS sources, the agreement between model predictions 
and observational data appears to be satisfactory. 

Substantially different evolution properties are re- 
quired for GPS galaxy and quasar populations, while the 
derived local luminosity functions of the two populations 
have a similarly flat slope. This similarity may, however, 
be fortuitous, particularly if the radio emission of GPS 
quasars (but not of galaxies) is boosted, and the lumi- 
nosity function flattened (Urry & Shafer 1984), by a rela- 
tivistic Doppler effect, as suggested by some observational 
data (O'Dea 1998). 

The model predicts that, if the peak frequency distri- 
bution extends to very large values (~ 1000 GHz in the 
rest frame), GPS sources (mostly quasars) may comprise 
quite a significant fraction (30 to 70%) of bright (S > 1 Jy) 
radio sources at v > 30 GHz, according to estimates by 
Toffolatti et al. (1998) and Sokasian et al. (1999). If the 
maximum intrinsic peak frequency is ~ 200 GHz [an in- 
trinsic v p > 190 GHz has been reported by Grainge & 
Edge (1998)] the number of GPS quasars brighter than 1 
Jy at 30 GHz decreases by about a factor of 4 and the 
total areal density of GPS sources at the same flux limit 
decreases by a factor of about 2. 

In any case, the counts of GPS sources are expected 
to be substantially flatter than those of other classes of 
extragalactic sources. As shown by Figs. 3 and 4, the 
slopes of integral counts of both galaxies and quasars 
never exceed 1 and are frequently substantially smaller. 
This result may be surprising in the case of quasars, given 
the strong luminosity evolution of their population. One 
should remember, however, that as z increases, the ob- 
servational selection picks up higher and higher values of 
intrinsic peak frequencies and earlier and earlier phases, 
when the luminosity is higher, of the evolution of indi- 
vidual sources [cf. Eq. (^)J. The duration of early phases, 
however, rapidly decreases with increasing luminosity; the 
effect is, in some sense, equivalent to a negative density 
evolution which compensates for the positive luminosity 
evolution of the quasar population. This "negative evolu- 
tion" effect is much weaker in the case of galaxies because 
their redshift distribution is weighted towards lower val- 
ues and because the evolution of individual sources is less 
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rapid (p = — r/ + aX ~ —4.2 for galaxies while p ~ —6.7 
for quasars). 

The flatness of counts obviously implies that the frac- 
tion of GPS sources quickly decreases with decreasing flux; 
as a consequence, their contribution to small scale fluctu- 
ations in the frequency range covered by Planck Surveyor 
and MAP missions is expected to be minor. 

In closing, we stress once again the exploratory nature 
of this study. Several conclusions are still tentative and will 
need to be reassessed as more data will become available. 
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